Designing reservoirs for 1/t decoherence of a qubit 
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Anomalous decoherence in the Jaynes-Cummings model emerges for a certain class of bosonic 
reservoirs, described by spectral densities with a band edge frequency coinciding with the qubit 
transition frequency. The special reservoirs are piecewise similar to those usually adopted in Quan- 
tum Optics, i.e., sub-ohmic at low frequencies and inverse power laws at high frequencies. The 
exact dynamics of the qubit is described analytically through Fox //"-functions. Over estimated 
long time scales, decoherence results in inverse power laws with powers decreasing continuously to 
unity, according to the particular choice of the special reservoir. The engineering reservoir approach 
is a new way of strongly delaying the decoherence process with possible applications to Quantum 
Technologies, due to the simple form of the designed reservoirs. 

PACS numbers: 03.65.Yz,03.65.Ta,03.65.-w 



I. INTRODUCTION 

Decoherence indicates the process that a quantum sys- 
tem of interest undergoes through the interaction with 
its external environment. The corresponding time evolu- 
tion and the destructive effects on quantum coherence are 
treated in the Theory of Open Quantum Systems [H [|| . 

Great attention has been devoted to the dissipative ef- 
fects of a two-level system (TLS), a qubit in Quantum 
Information Theory, interacting with an external envi- 
ronment modeled by a reservoir of bosons The 
applications of this simple model are most various: from 
Nanotechnology to Quantum Information and Quantum 
Computing, from Quantum Optics to circuit QED, to 
name a few. Still, the central issue and one of the great- 
est challenges remains the way to control or delay the de- 
structive effect of the external environment on coherence. 
For example, the decoherence time @ in Magnetic Reso- 
nance has orders of magnitude ranging between nanosec- 
onds and seconds. 

Various techniques are adopted in order to give an an- 
alytical description of the exact dynamics of the qubit. 
Interesting results emerge from the adoption of the re- 
solvent operator Q in rotating wave approximation 
with a Lorentzian distribution of field modes The 
assumption that the coupling constants vary slowly with 
frequency, allows a complete analytical treatment and the 
exact dynamics results in oscillating behaviors enveloped 
in exponential decays. For a detailed report we also refer 
to [9]. 

An interesting model for the spontaneous decay of a 
TLS in a structured reservoir, has been introduced by 
Garraway [H| and solved exactly for a generic dis- 
crete or a continuous distribution of field modes, de- 
scribed by Lorentzian type spectral densities and a spe- 
cial non-Lorentzian one with two poles in the lower half 
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plane. The crucial conjecture that the frequency range is 
(— oo, +oo), allows the analytical evaluation of the exact 
dynamics in terms of the pseudomodes [l2| and the poles 
of the spectral density in the lower half plane. The time 
evolution results in oscillations enveloped in exponential 
relaxations. Recently [la ], the model adopted by Gar- 
raway has been chosen to compare the qubit dynamics of 
an exact master equation in time convolution less form 
with the Nakajima-Zwanzig master equation, through a 
perturbation expansion of the memory kernel. 

The realization of structured environments providing a 
discontinuity in the distribution of the frequency modes, 
also named photonic band gap (PBG) [IJ-QtI; introduces 
new phenomena in the atom-cavity interactions. For ex- 
ample, the spontaneous emission of a two-level atom near 
the edge of a PBG exhibits oscillatory relaxations [18j in- 
stead of a purely exponential decay. A theoretical model 
providing a PBG structure in a A^-period one dimensional 
lattice has been proposed in Ref. [l9| by arranging an 
appropriate sequence of the unit lattice cells. The den- 
sity of the frequency modes is analytically evaluated as a 
function of the transmission coefficient of each unit cell. 

In line with the attempt to contain the destructive ef- 
fect of the external environment on a qubit, decoherence, 
a special reservoir of bosons is designed in Ref. [20| with 
a PBG structure, excluding modes with frequencies lower 
than the transition frequency of the qubit. The exact dy- 
namics is described analytically by a linear combination 
of incomplete Gamma functions and decoherence results 
in a 3/2 inverse power law relaxation over an evaluated 
long time scale. If compared to the exponential like re- 
laxations, decoherence is strongly delayed. In this sce- 
nario, we adopt the Jaynes-Cummings model, starting 
from the initially unentangled states of the system qubit 
plus reservoir, adopted in Ref. [To| . and design a wider 
class of reservoirs inducing slower forms of relaxation for 
a further delay the decoherence process. We notice in 
advance that the analytical calculations leading to the 
exact dynamics, rest on a continuous distribution of pos- 
itive mode frequencies. 
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II. THE MODEL 



The interaction between the qubit and the dissipative 
environment is described through the Jaynes-Cumming 
model with a continuous distributions of field modes [1, 
0, IE E3, EH ■ By choosing H = 1 , the Hamiltonian of the 
whole system is Hg + He + Hj , where 

oo 

H s = wocr+o--, H E = y^ oj k q[q fc , 



fc=i 



g k cr + <K> a fc + 3 fc cr_ <K> a*. 



fe=i 



The rising and lowering operators, ct+ and <r_, respec- 
tively, act on the Hilbert space of the qubit, defined 
through the equalities a+ = <rl_ = |1)(0|, while a\ and 
a k are the creation and annihilation operators, respec- 
tively, acting on the Hilbert space of the /c-th boson, 



fulfilling the commutation rule ak,a' k = S k ^ k i for ev- 
ery k, k' — 1,2,3,.... The constants g k represent the 
coupling between the transition |0) <H> |1) and the k-th 
mode of the radiation field, while ujq is the qubit transi- 
tion frequency. In the following we refer to the system 
of a TLS interacting with a cavity supplying a reservoir 
of field modes, as studied by Garraway [To[ and adopted 
in Ref. [l3j]. Starting from the initial state of the total 
system 



|*(0)) = (co|0)+ci(0)|l))®|0> B , 



(1) 



where \0)e is the vacuum state of the environment, the 
exact time evolution is described by the form 

|*(i)> = c |0) ® |0} B + ci(i)|l) ® |0) B 

OO 

+ ^6 fc (t)|0>® \k) E , |fc) E = 4|0) E , fc = 0,l,2,.... 
fc=i 

The dynamics is easily studied in the interaction picture, 
|*(t))j - e*< Hs+fl *)*|tf(t)) = c |0) ® \0) E 

OO 

+ Ci(t)|l) <8> |0> B + X)Bfc(t)|0) ® 
fc=i 

where i is the imaginary unity, C\(t) = e* Wot ci(t) 
and B k {t) = e lulkt b k (t) for every k = 1,2,.... The 
Schrodinger equation gives the forms: 



Cl{t) = -iJ29kB k {t), 



-%(u) k —u )t 



fc=l 



leading to the following convoluted structure equation for 
the amplitude (1| <g> B (0||*(t))/, labeled as Cx(t), 



where / is the two-point correlation function of the reser- 
voir of field modes, 



f(t-t')=J2\9k\ 2 e- l{ ^ o) ^ t 'y 
fc=i 

For a continuous distribution of modes described by 
r)(uj), the correlation function is expressed through the 
spectral density function J (to), 



/(r)= / J(oj)e-<"-"°>duj, 
Jo 

where J (uj) = r] (uj) \g (uj)\ 2 and g(uj) is the frequency 
dependent coupling constant. 

The exact dynamics of the qubit is described by the 
time evolution of the reduced density matrix obtained 
by tracing over the Hilbert space of the bosons, 

Pi,i(*) = l-fl>,o(*)=Pi,i(0) |G(i)| 2 , (3) 
Pi,oW=pS > i(*)=Pi,o(0)e- Wo *G(t). (4) 

The function G(t), fulfilling the convolution equation 

G(t) = - (/ * G) (t), G(0) = 1. (5) 

The function drives both the dynamics of the levels pop- 
ulations and the decoherence term. 



III. THE EXACT DYNAMICS 

We study the exact dynamics of the reduced density 
matrix of the qubit, interacting in rotating wave approx- 
imation with a reservoir of bosons described by the con- 
tinuous spectral density 

T i s 2A (uj - uj a ) a 6 (uj - UJp) 
Ja{UJ) = ; -J , (6) 

a £ + (uj - ojq) 
A > 0, a > 0, 1 > a > 0. 

This simple form exhibits a PBG edge in the qubit tran- 
sition frequency, has an absolute maximum M a at the 
frequency Q a , 

M a = J a (tt a ) = Aa a ' 2 a a - 2 (2 - a) 1 ' 01 ' 2 , 
fl a = uj + aa 1/2 (2-a) 1/2 . 

The above spectral densities are piece-wise similar to 
those usually adopted, i.e. sub-ohmic at low frequen- 
cies, uj ~ wo, and inverse power laws at high frequencies, 
uj 3> ujq, similar to the Lorentzian one, though with dif- 
ferent power, 



C 1 (t) = -(f*C 1 )(t), 



(2) 



J a (uj) ~ 2 Aj a 2 (uj — ujo)° , 
J a (uj) ~2Auj a ~ 2 1 uj^ 



-co. 
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The exact dynamics of a qubit interacting with a reser- 
voir of bosons described by the spectral density J a (w), 
is driven by the function G a (t) , solution of Eq. (JSJ , 



{-\) n n\zl z™~ k t 3n - ak 



n=0 k=Q 



k\{n-k)\ 

x (E^ n _ ak+1 (~ Zl t 2 ) - a^E^_ ak+3 (- Zl t 2 ) ) , (7) 

expressed as a series of Generalized Mittag-Leffler func- 
tions SI El, 



a, ft 7 G C. 



^ n!r(an + /3)' 

n— 



9?{a}>0, 3?{/3}>0, 



where (7) = 1 and (7) = T (7 + n) /T (7). The param- 
eters involved read 

Z\ = -KAa a ' 1 sec (na/2) — a 2 , zq = inAa a cos (7ra/2) , 



z Q = -2z7rAe" I7rQ/2 csc(7ra) . 



(8) 



The proof is performed below. 

The Generalized Mittag-Leffler function, fundamental 
in Fractional Calculus [2lj . is a particular case of the 
Fox .ff-function, defined through a Mcllin-Barnes type 
integral in the complex domain, 



Trm.n 



(ai, ai) , . . . , (a p , a p ) 
(6 1 ,/3 1 ),...,(6 9 ,/3 9 ) 



n^r (6j + gys) n^r (1 - ai - z~ 6 
nf =n+1 r (a, + a iS ) n? =m+1 r (1 - b 3 - p jS \ 



1 

2th 
■ ds. 



Under the conditions that the poles of the Gamma func- 
tions in the dominator, do not coincide, also the empty 
products are interpreted as unity. The natural numbers 
m,n,p,q fulfill the constraints: < m < q, < n < p, 
and a>i,/3j 6 (0, +00) for every i — 1, ••• ,p and j = 
1, ■ • ■ ,q. For the sake of shortness, we refer to [23j for 
details on the contour path C, the existence and the prop- 
erties of the Fox ii- functions. The relation 



1 



r(7) 



#1,2 



(1-7,1) 
(0,1), (1-/3, a) 



leads to a series solution of Fox H functions 

(-i) B *5*o 



k ^n — k j.3n — ak 



ri=0 fc=0 
2 



fc!(n-fe)! 

(-n,l) 
(0,l),(afc-3n,2) 



(-n,l) 

(0, 1) , (afc - 3ra - 2, 2) 



(9) 



The Generalized Hypergeometric, the Wright [24[ and 
the Meijer G-functions [? ] are particular cases of the 



Fox if- function, thus, G a (t) can be expressed as a series 
of each of these Special functions, as well. 

Particular cases give simplified solutions. For example, 
the condition A = A w , 

,4 W = cos (wa/2) , (10) 

7T 

corresponding to Z\ = 0, gives a power series solution, 

v n I 1 \n „| „fc ~n— j.3n-afe 

« u jfc? ( n _ r ( 3n _ afc + 1) 

\ , r(3n-afc+l) ,1 
X { 1 -° r(3n-afc + 3) * }' 1 > ° > °" ™ 

If the parameter a takes rational values, p/q, where p 
and q are distinct prime numbers such that < p < q, 
the solution of Eq. ^ can be expressed as a modulation 
of exponential relaxations [26| . 

/>oo />oo 

<W*) = / rf77 / d£<W^0e-«', (12) 
Jo Jo 



where 



n nil , s 



=1 fc=i 



x sin 1 



The rational functions 6;^ (z) read 

d k ~ x ldz k - x [(z« - a) (z" + a)(z- Ci) m ' IQ (z)] 



h,k (z) 



(mi - k)\ (k - 1)! 



for every I = 1, ■ • • , n, and fc = 1, . . . , m;. The complex 
numbers Cw " , Cn are t ne roots of the polynomial 



Q(z) = Z 3q + Z lZ q + Z a Z P + Z Q 



(13) 



and mi is the multiplicity of Q, for every I — 1, . . . ,n, 
which means Q(z) = n™ =1 (z — Q) m ' and m i = 3q- 

The case a = 1/2 exhibits a simplified exact dynamics 
described by a finite sum of Eulerian functions [20j . For 
a = 3/4 and A = a 9 / 4 cos (37r/8) /71", the parameter Z\ 
vanishes and the roots Ci 1 ■ ■ • 1 can b e evaluated ana- 
lytically from the solutions of a quartic equation. We do 
not report the expressions for the sake of shortness. In 
the remaining cases of rational values of a, the roots of 
Q(z) must be evaluated numerically, once the numerical 
values of both A and a are fixed. These details complete 
the necessary analysis of the function G a (t). 

Finally, the exact time evolution of the qubit is ob- 
tained from Eqs. ([3]) and ((3]), by replacing the function 
G(t) with G a (t), analyzed in the present Section. 
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IV. INVERSE POWER LAWS 

The theoretical analysis of the exact dynamics, per- 
formed above, leads to the following concrete result: a 
time scale r emerges such that, for t ^> r, the function 
G a (t) exhibits inverse power law behavior described by 
the asymptotic form 

G a (t) Vat- 1 - 01 , i^+oo, l>a>0, (14) 

where 

2ma 2 ( 1 - Q )e- ,TO / 2 esc (ira) sec 2 (na/2) 



V a = 
A simple choice is 
t = max < 1 , 



7rAr(l - a) 



1/3 



1/a 



(15) 



the proof is performed below. 

Thus, over long timescales, t 3> t, the qubit exact 
dynamics is described by inverse power law relaxations: 

p M (t) = 1 - Po ,o(t) ~ Pi,i(0) \V a \ 2 t- 2 - 2a , (16) 

(17) 



Pi,o(*) = PoAt) = pi,o(0)P Q e- lWot t 



— tujQt j.—1—a 



for every a 6 (0, 1). 

In Ref. the decoherence processes corresponding 
to a Lorentzian and J\i% (uj) spectral densities are com- 
pared. The exponential type relaxations emerging in the 
Lorentzian case, vanish faster than the inverse power 
laws related to the spectral densities J a (uj), for every 
a 6(0,1). 

The results obtained are summarized as follows. Start- 
ing from the initial condition (TTJ where the reservoir, in 
the vacuum states, and the qubit are unentangled, we 
evaluate analytically the exact dynamics of the qubit in- 
teracting, in a rotating wave approximation with a reser- 
voir of bosons described by the spectral density ©. The 
time evolution is expressed by the series ([9]) of Fox H- 
functions, and, over long timescale, {> r, decoherence 
results in an inverse power law relaxation proportional to 
f or ever y a £ (o, 1), according to the choice of the 
special reservoir ([6]). The qubit ultimately collapses into 
the ground state. 



V. CONCLUSIONS 

Anomalous forms of qubit decoherence emerge from 
the Jaynes-Cummings model for reservoirs of bosons de- 
scribed by special continuous spectral densities with a 
PBG edge coinciding with the qubit transition frequency. 
The designed spectral densities are piece-wise similar to 
those usually adopted, i.e., sub-ohmic at low frequencies 
and inverse power laws at high frequencies, similar to 
the Lorentzian one. Initially, the qubit and the reservoir, 



versing in the vacuum state, are unentangled. The ex- 
act dynamics is described analytically through series of 
Fox iJ-functions. Over estimated long time scales, qubit 
decoherence results in inverse power law relaxations with 
powers decreasing continuously to unity, according to the 
choice of the special reservoir. 

An environment supplying the designed reservoir of 
field modes can in principle be realized with materials 
providing the PBG structure. For example, a ./V-period 
one dimensional lattice can reproduce a band gap by ar- 
ranging the appropriate sequence of dielectric unit cells. 
The corresponding density of frequency modes is struc- 
tured by the transmission coefficients of each unit cell. 
Also, the advanced technologies concerning diffractive 
grating and photonic crystals allows the realization of 
tunable ID PBG microcavities [27l . [28| . The simple form 
of the designed reservoir may be accessible experimen- 
tally. The action of such a structured environment on a 
qubit could be a way of delaying the decoherence process 
with fundamental applications to Quantum Information 
processing Technologies. 



VI. PROOF 

A detailed demonstration of the solutions driving the 
exact dynamics follows. The reservoirs of bosons de- 
scribed by the following class of non negative, non di- 
vergent and summable spectral densities are considered: 



J (uj) du < oo, J (u>) = Q (uj — ujq) A (uj — luq) 



in this way, Eq. ([5]) gives 



G(u) = {u-iS(A)(-iu)y 1 . 
3?{u}>0, |arg{-zu}| < n, 



(18) 
(19) 



where S is the Stieltjes transform. The constraints: 
Jq 00 A (ui) du < oo, and 3? {u} > 0, guarantee the uniform 
convergence of integrals involved in the Integral trans- 
forms, so that the above equality holds true. The class 
of spectral densities ^ and Eq. §5§ lead to the following 
Laplace transform: 



G a (u) 



Z\U + Z a U L 



(20) 



the parameters Zq, z a and z\ are defined by Eq. ([8]). 

The function G a (t) is obtained through the convergent 
term by term Laplace inversion of the series expansion of 
Eq. (2D), 



- Z\U + z a u 
ik-n-l( a ,2 



(u 2 - a 2 ) 



)!=0 fc=0 

z a u a + z 



(u 2 + zi) 



n+l 



U A + ZiU 



k\(n-k)\ 

< 1. 
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This way, a series solution of Eq. ([5]) can be built in 
terms of either Generalized Mittag-LcfHer functions, Eq. 
(O , or Fox H functions, Eq. (j9|), for details we refer to 

mil. 

The inverse power law behavior of G a (t) over long time 
scales, is found through the series expansion 



oo n k 

EEE-r 



(-l) n n! 



u 3 + zm + z a u a + z ^f^f^jKn ~ k )K n - i) ! 



- n -l z n-k z k- j{u 2 



,2\ .Mn-k)+k+2j 



holding true under the constraint 



u + Z\u + z a u c 



< 1. 



(21) 



(22) 



The formal term by term Laplace inversion of Eq. (|2"Tj) 

gives 



oo n k 



n=l fe=0 3=0 



71-1 



-fe 



x a 



j\ (n - k)\ (n - j) ! 
2 T (a(fc — n)—k — 2j — 2) _ 



T(a(k -n)-k- 2j) 

j.—l—(x{Ti—k)—k—2j 

-, t^+oo, (23) 



r(a(fc-n)-fc-2j-2) : 



leading to the asymptotic solution (1141) . The time scale 
for inverse power law behavior descends from the inequal- 
ity (f22j) and the constraint t ^> 1, requested by Eq. ([23]). 
A possible choice is obtained by imposing that the abso- 
lute value of each term of the left side of the inequality 
([22]) is less than 1/3. This way, an upper bound for \u\ is 
obtained and the corresponding inverse estimates a time 
scale for inverse power laws, given by Eq. (|15[) . 
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